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I. INTRODUCTION 

The statistical characterization and modeling of two- 
phase disordered microstructure is a central problem in 
many fundamental and applied sciences jjj. Predicting 
the properties of disordered materials relies on the avail- 
ability of accurate microstructural models, which rely in 
turn on accurate statistical characterization. After the 
volume fraction of each phase, and interfacial surface 
area, the most important statistical quantity is the two- 
point correlation function which is obtained from cross- 
sectional micrographs, or small angle scattering exper- 
iments. Although the two-point correlation function is 
very useful, there are a variety of important problems 
where more detailed statistical information is necessary. 
Another useful characteristic of microstructure, which 
has proved essential in theory and application, is the 
chord-length distribution function 0]. The chord func- 
tions play an important role in stereology ||, mineral- 
ogy B, the interpretation of small angle x-ray scattering 
data Tq], and have been incorporated in theories of mass 
transport in porous media Recently the chord func- 
tions (and the related "lineal-path" function 0) have 
been employed in the generation of three-dimensional 
(3D) microstructural models for predicting macroscopic 



properties |^-[Io|. In this paper, we derive approximate 
forms of the chord functions for a relatively new model 
of random media based on Gaussian random fields. 

A useful model of two-phase random porous and com- 
posite media is obtained by modeling the internal inter- 
face of the microstructure as the iso-surface (or level-cut) 
of a correlated Gaussian random field y(r) |Tl]-[lq]. A re- 
gion of space can be divided into two phases (e.g. pore 
and solid) according to whether y(r) is less or greater 
than some threshold. We define phase 1 to occupy the 
region where y(r) < j3 and phase 2 to occupy the region 
where y(r) > [3. This is illustrated for two dimensions in 
Figs. l](a) and (b). The model, and variants, have proved 
useful in describing the microstructure of many different 
materials j^|l(||l7|] ; and a more thorough characteriza- 
tion of the microstructure is an important goal. The 
two- and three-point correlation functions of the model 
can be calculated, but the chord functions are presently 
measured from simulations [00] • There are advantages 
(for speed, accuracy and interpretation) in obtaining ana- 
lytic expressions for the chord functions of the model. As 
we show, the problem is equivalent to finding the proba- 
bility density governing the time it takes for a Gaussian 
random process to first cross an arbitrary threshold. This 
is a conventional first-passage time problem. 
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FIG. 1. Generation of a two-phase model (b) by threshold- 
ing a Gaussian random field (a). The chord length distribu- 
tion functions are calculated by counting the number of chords 
of a given length (b). In one dimension the chord lengths z\ 
and 22 are defined by up- and down-crossings (shown as u and 
d respectively) of a threshold j3 by a random process (c). 



First-passage time problems arise in many branches of 
physics |XS| , information theory ||{J , queuing theory , 
ocean science jElJ and reliability studies in the engineer- 
ing sciences [|22| amongst others. For this reason the 
problem has received a great deal of attention. Rice 
actually provided a formal series solution to the prob- 
lem for Gaussian processes ^3j. However, the series 
involves very difficult integrals of which only the first 
is generally evaluated |i~8[ ; the results being accurate 
for small time. There are many approaches to finding 
useful approximations for first-passage times (the afore- 
mentioned references provide reviews of the literature 
in each field). Here we restrict attention to methods 
based on the assumption that the lengths of successive 
chords are independent. This idea can be traced back to 
Siegert |24| and McFadden [ p5[ . The approximation we 
use is most clearly derived from the independent interval 
process |5|,p6fl for which the assumption of independent 
chords is true by definition. This approach was recently 
suggested for the calculation of persistence times in coars- 
ening p7lp8| ■ The method can be extended to obtain the 



chord-distributions, and actually provides an extremely 
useful way of viewing related problems in different fields. 
In the following sections we derive some properties of the 
chord function and the independent interval process. The 
model is then applied to approximate the chord functions 
of level-cut Gaussian random fields. Finally, we compare 
the theoretical results with experimental data. 



II. CHORD DISTRIBUTION FUNCTIONS 

For a two-phase medium, there is a chord-distribution 
associated with each phase Pi(z) (z=l,2). The quantity 
Pi(z)dz is defined as the probability that a randomly cho- 
sen chord in phase i (a line segment with end points on 
the phase interface) has length in the range [z, z + dz]. In 
Fig. 0(b), we illustrate calculation of pi(z) andp2(z) for 
isotropic media: An "infinite" line (or many shorter ones 
with arbitrary orientation) is drawn through the mate- 
rial, and the number of chords in phase 1 [iVi(z)] with 
length in the range [z, z + dz] is counted. If N is the total 
number of chords of phase 1, then pi(z)dz — N±(z)/N. 
The quantity P2 is defined in an analogous manner. An 
important quantity is the number of phase crossings per 
unit length n c = 2N/L. A fundamental relation in stere- 
ology gives n c — ^s v where s v is the specific surface, 
i.e., the surface area to total volume ratio of a 3D com- 
posite. Another useful relation follows from the fact 
that total length of the chords is equal to that of the 
original line: Yl z Ni{z)z + ]V 2 (z)z = L. If we di- 
vide by N and convert the sums to integrals we have 
/ p\(z)zdz + J pi(z)zdz = 2/n c . Similarly ^ 2 Ni(z)z = 
faL where fa is the volume fraction of phase i. This 
gives the relation (zi) = J pi(z)zdz — 2fa/n c . These 
equalities are used extensively below. A useful statistical 
quantity is the lineal-path function of each phase Li{z) 
(i = 1,2), which represents the probability that a ran- 
dom line segment of length z thrown into the material 
falls completely within phase i and is related to Pi(z) by 
Li(z) = \n c J z °° Pi(x)(x - z)dx 0. 

One reason for the usefulness of the chord function is 
that pi (z) can be directly interpreted in terms of observ- 
able microstructure features. First, if pi(0) ^ then, at 
the resolution at which it is measured, phase i contains 
sharp corners. Second, the value at which Pi(z) takes a 
maximum value provides an estimate of the length scale 
associated with phase i. Another is provided by the aver- 
age chord length (zt). Third, if pt(z) ^ for large z, con- 
nected regions in phase i at scale z exist. This "connect- 
edness" information (along a lineal path) is clearly im- 
portant if long-range phenomena (like macroscopic prop- 
erties) are to be studied. The direct relationship between 
Pi (z) and morphology shows that the chord functions give 
a strong "signature" of microstructure and are therefore 
an important tool in the characterization of porous and 
composite media. 

The chord and lineal-path functions are closely related 
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to first-passage times in the theory of stationary time- 
dependent random processes. The analogy is shown in 
Fig. 0(c). If the values of a random-field along a line 
are plotted against time a one-dimensional (ID) random 
process y(t) is obtained. If a down-crossing of the thresh- 
old (3 occurs at t = td [so y(td) — (3), the probability 
that the process first exceeds the threshold in the interval 
[td + t, td + 1 + dt] is pi{t)dt. The corresponding density 
for the first down-crossing after an up-crossing (some- 
times called the second-passage time) is P2 (t) . There are 
several other common first-passage times which relate to 
our work. Suppose y(t) is a random process represent- 
ing the response of an electrical or mechanical component 
which fails iiy(t) exceeds some (generally high) threshold 
(y(t) > (3). A key quantity is the probability of failure 
F(t) in the interval [0,i] [E9[. Since the safe region is 
phase 1 the failure probability is just F(t) = 1 — Li(t) 
[F(0) = 1 — Li(0) = 02 being the probability of instan- 
taneous failure]. A related quantity is the probability 
distribution of failure times fi(t) given that t = is in 
the safe region (22). This is given by fi{t) = — L[(t) / . 



III. INDEPENDENT INTERVAL PROCESS 



The independent interval process 26 is constructed 
by taking the lengths of successive intervals to be inde- 
pendent of one another and distributed according to the 
probability (chord) distribution functions p±( z) andp2{z) 
for phase 1 and 2 respectively (this implies a stationary 
process). Since this completely defines the process it is 
possible to derive all other statistical properties (such as 
the two-point correlation function) frompi(z) andp2(z). 
The indicator function /(x) (which is unity in phase 1 
and zero in phase 2) is very useful in this regard. From 
above we have 



= (inx)i) = 



(J(x)) - -n c ( Zl ) 

(1 - Ax)) = ~n c (z 2 >. 



(1) 
(2) 
(3) 



Note that the results are true irrespective of correlations 
between successive chords. 

Now consider the two-point correlation function which 
is defined as 



S u (x) = (/( Xl )/(x 2 )). 



(4) 



Stationarity of the process (and isotropy in more than 
one dimension) implies that the average only depends on 
x = |x2 — xx|. To express Su as a function of p± and P2 
it is necessary to derive some preliminary results. The 
method is adapted from Ref. p7[ . Although the model 
is independent of any random field it is useful to de- 
scribe the left-hand ends of chords in phase 1 and 2 as 



down-crossings and up-crossings, respectively (See Fig[l]). 
Consider the probability R u i(x) that a point at distance 
x to the right of an up-crossing (or right-hand end of a 
chord in phase 1) falls in phase 1. Suppose the chord (of 
phase 2) immediately to the right of the down-crossing 
has length z 2 . If Z2 > x then the point falls in phase 2. 
If Z2 < x the chance that the point falls in phase 1 is 
Rdi{x — Z2), where Rdi{y) is the probability that that a 
point a distance y to the right of a down-crossing falls in 
phase 1. Thus 



R u i{x) = / dz 2 p2 (z 2 )Rdi (x - z 2 ). (5) 
Jo 

Now consider the converse problem for Rd\(x). Let z\ be 
the length of the first chord (of phase 1) to the right of a 
down-crossing. The probability that a point a distance x 
from the down-crossing falls in phase 1 is unity if z\ > x 
and R u i{x — z\) if z\ < x giving 



dzi pi(zi) + / dzi pi(zi)R ul (x - zi). (6) 



Rdi (x) 



Taking Laplace transforms of both equations and solving 
we have 



Rui 

Rdl 



P2O- -Pi) 
s(l - P1P2) 

s(l -P1P2) 



(7) 
(8) 



where f=f(s)= J °° e~ sx f(x)dx denotes the usual Laplace 
transform. 

Next we need the probability Qi(y)Sy that an arbi- 
trary point will fall in phase 1 a distance [y, y + Sy] from 
the left of the first down-crossing on its right. The chance 
that the point falls on a 1-chord of length [z\, z\ + dz{\ is 
\n c p\(z\)z\dz\. The point will be uniformly distributed 
along the chord, therefore if it falls on a chord of length 
Zi the probability that it lies a distance [y, y + Sy] from 
the left end is just 5y/z\. To obtain the total probability 
we must sum over all chords with z > y giving, 



n c f°° Sy 

2 Jy z\ 

Q^d-Pr). 



(9) 
(10) 



Note that Qi(0) = \n c as it should since the probability 
that a point lies within a distance Sy (to the left) of a 
down-crossing must be ^n c Sy. 

Recall that Si 1(2) is the probability that two points 
a distance x apart fall in phase 1. The chance that the 
second point falls in phase 1 (given that the first does) 
depends on whether the distance to the first up-crossing 
y is greater or less than the distance x. If y > x then the 
right-hand point falls in phase 1 with unit probability. 
If y < x then the right-hand point falls in phase 1 with 
probability R u i(x — y). Thus we have 
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Su(x) 



dyQi{y)+ / dyQ 1 {y)R ul {x-y). (11) 
Jo 



Note that <Sii(0) = <fii as it should since 5 l n(0) is just the 
probability that a single point falls in phase 1. Taking 
Laplace transforms we obtain 



a _ 4>i n c (1 -pi)(l -p 2 ) 

OH — 7T2 i ~ ' 

s 2s z 1 - pip 2 



(12) 



Two other correlation functions of the independent in- 
terval process are important for later discussions. The 
first, S c i(x) = (|/'(xi)||/(x 2 )|) is the ID analog of the 
surface- void correlation function which arises in the study 
of 3D porous materials 0. The second is the crossing- 
crossing correlation function S cc (x) — (|/'(xi)| |/'(x 2 ) ) 
analogous to the surface-surface correlation function ^_\. 
Again, stationarity implies that S c i and S cc depend only 
on the distance x = |x 2 — xi|. 

From the definition of S c i we have 



eS c i (x) 



J(xi 



/(xi 



(13) 



where e is small. The expression on the right-hand side 
is the probability that X\ lies within a distance | (which 
we call an e-interval) of a crossing and that x 2 lies in 
phase 1. Without loss of generality, we consider the case 
X2 > xi. Now there is an equal chance of the first point 
landing in an e-interval of an up or down crossing (prob- 
ability ^en c ). The probability that x 2 falls in phase 1 is 
then either R u \(x) or R<n{x). Hence, for e — > 0, 



eS c i(x) = ien c x [R u i(x) + R d i(x) 
or, using Laplace transforms, 



Sci — 



n e (l-p 1 )(l+p 2 ) 
2s 1 - pip 2 



(14) 



(15) 



Similarly we write the crossing-crossing correlation 
function as 



e 2 S cc (x) 



I(*2 + -e) - Z(x a - -e] 



(16) 



where the expression on the right is clearly the probabil- 
ity that both points lie within an e-interval of a cross- 
ing. To express this quantity in terms of the chord- 
distributions several additional functions are needed. Let 
Rud{x) be the probability that a point a distance x from 
an up-crossing falls in an e interval of a down crossing. 
This can occur in two ways; either the first chord of phase 
2 adjacent the up-crossing has length z 2 = x, or z 2 < x 
in which case the point falls near a down crossing with 
probability Rdd(x — z 2 ), where Rdd{y) is the probability 
that a point a distance y to the right of a down-crossing 



is in an e-interval of a down-crossing. This is expressed 
as 

Rud{x) — ep 2 (x) + / dz 2 p 2 (z 2 )R d d(x ~ z 2 ) (17) 
Jo 

The functions R uu , Rdu are similarly defined and three 
additional relations amongst the four functions can be 
derived and solved to give 

Rud = efc/(l-PiP2) (18) 
Rdu = ep 1 /(l-pip 2 ) (19) 
Rdd = Ruu = epifc/ (1 - pipi)- (20) 

Now from the definition of S cc it is clear that 



e 2 S cc {x) = -en c [R uu + R ud 



Rdu + Rdd] (21) 



as e — > 0. After taking Laplace Transforms this gives 



S c 



n c P1+P2+ 2fflip 2 
2 1 - pip 2 



(22) 



A final expression which is useful for relating this work 
to prior approaches is 

Ri u (x) = efi(x) + / dyR uu (x - y)fi(y), (23) 
Jo 

where Ri u (x) is the probability that a point a distance 
x from a point in phase 1 falls within an e-interval of an 
up-crossing, and /i(x) is the probability distribution of 
"failure times" discussed earlier. 

Although this analysis of the independent interval pro- 
cess, and the results for S c i and S cc appear new, variants 
of Eq. ([TJ|) for Su have actually been derived in several 
quite different contexts. It is useful to briefly demon- 
strate these connections. The Poisson-Boolean model 0] 
of random two-phase media has been widely studied and 
applied. In this model, grains (which may be of differ- 
ent shapes and sizes) are placed at uncorrelated random 
points in space (so the grains may overlap) . For spherical 
inclusions this it is just the overlapping sphere (or Swiss- 
cheese) model [0. For convex grains of any shape or 
size, the chord-distribution function of the phase exterior 
to the grains is Pi(z) = Aexp(— Xz) where A = s v /A<f>i. 
Since the grains are uncorrelated in space and convex it 
is clear that the length of successive chords along any line 
through the a realization of the model will be indepen- 
dent. Therefore Eq. ( p"2| ) applies. Specializing to Boolean 
models [pi = A/(A + sjj we recover the well known result 
(e.g. Ref. |) 



sS u 



1-Pi 
1 -P1P2 



(24) 



Since Sii(x) is known for many Boolean models, this re- 
sult allows p 2 to be calculated; 



p 2 = 1 + s/A- 4>i/XSu. 



(25) 
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Let us further specialize to the case of a ID Boolean 
process where the grains are rods with random 
lengths distributed according to the cumulative dis- 
tribution function ^>(z). For this model Sn(x) — 
4>i exp (—A Jq(1 — ^(z))dz) and the formula for p 2 be- 
comes a well known result for the busy period in an 
M/G/oo queue 

The independent interval process is useful in the in- 
terpretation of small-angle x-ray scattering data. In 
Refs. f||2(| the Fourier transform of d 2 Sn/dx 2 was de- 
rived in terms of the Fourier transforms of the chord 
functions. The result can be shown to be equivalent to 
Eq. (|l2|). Rice's j2^] formulae for Sn for two different 
types of random telegraph signal can also be re-derived 
using Eq. (|l2|). In the next section we show how the 
independent interval process can be used to derive use- 
ful approximations of the chord functions for the single 
level-cut GRF model. 



and their spatial derivatives (Vy(rj), etc.) are correlated 
Gaussian random variables with known joint probability 
distribution. This allows many useful statistical proper- 
ties of the thresholded model (such as Sn) to be calcu- 
lated. 

In the previous section we derived the properties of 
the independent interval process in terms of p\ and p 2 . 
However, it is clear that if any two of the microstructure 
functions (pi, p 2 , Sn, S c i, S cc ) are known the remainder 
can be determined by means of Eqs. (p"2|), ( |l5| ) and (p^). 
Since Sn, S c i and S cc can be calculated for the level-cut 
GRF model, this allows an approximation to be derived 
for the chord functions. The accuracy of the results will 
then depend on the validity of the hypothesis that succes- 
sive chords of the model are uncorrelated (or nearly so). 
For simplicity we consider the case where Sn and S c i 
are known. Simultaneous solutions of Eqs. (]l^) and ( |T5| ) 
then give 



IV. CHORD FUNCTIONS OF LEVEL-CUT 
GAUSSIAN RANDOM FIELDS 

A Gaussian field is statistically specified in terms of 
a field-field correlation function (y(r 1 )y(r 2 )) = g(|r 2 — 
ri | ) . Here we consider isotropic stationary fields 

with zero mean [(y(r)) = 0] and unit variance [(y 2 (r)) = 
g(Q) = 1]. Many methods exist for generating GRF's 
with a given g(r). For example, in ID we have 



N 

E« 

1=1 



2tti 

,; cos kiX + bi sin kix; fcj = (26) 



where di and bi are independent Gaussian random vari- 
ables with mean zero and (a 2 ) = (b 2 ) — ^ x F(ki). In 
the limit N, T — > oo such that N/T — > oo the correlation 
function is g(x) = F(k) cos(kx)dx. F{k) is called a 
spectral density. When F(k) is a wide or narrow distri- 
bution, the process is, respectively, said to be wide- or 
narrow-band. The 3D analog is 



JV 



N 



vw= E E E 



(27) 



l=-N m=-N n=—N 



where ki mn = ^{li + mj + nk) and Q, m ,„=a; mn + ify„ m . 
For y real and (y) — we take ci tm-n =C-i t - m ^ n and 
Co. o,o = 0. As above, ai mn and bi mn have mean zero and 

( a LnM & L«>=3 (¥) 3 x Pihmn)- In this case g(r) = 
Airr~ l J Q 4irkp(k) sin kr dk. A two-dimensional random 
field is shown in Fig. |l](a). 

A two-phase level-cut GRF model is specified by the 
microstructure indicator function I(r) = H(f3 — y(r)), 
where H is the Heaviside step function. There are two 
very useful properties of the Gaussian model. First the 
random field is ergodic (ensemble averages equal spatial 
averages), and second the variables y(ri) i = 1,2,3... 



Pi 



P2 



s(S c i - sSn 



n c - s(S c i + sSn - fa) 

Sci + 5-5*11 - 4>i 



S c i — son + ' 



(28) 
(29) 



This result can be considered a generalization of a recent 
approximation developed independently by the authors 
of Refs. [^?J and |28|. In a study of the zero-threshold 
case [where p{z)=pi(z)=p2(z)\ they found 



P = 



n c + s(2S n - 1) 
n c - s(25n - 1) 



(30) 



where Sii(x) = \ + ^ &rcs'm[g(x)]. This result is ob- 
tained by substituting ^i = 02 = | and S c i = n r ./2 
(which is true for any symmetric medium) into Eqs. (|28| ) 
and (|2^). The result provided an excellent approximation 
for the case g(x) = [1/ cosh(x/2)] d / 2 (d = 1,2,3). 

Clearly other approximations can be obtained using 
the independent interval process, and several have been 
previously given. For example Eq. ( |l7|) was obtained 
by McFadden |^5|] and Rice |3^], and approximate forms 
of Eq. (|23|) have been used to obtain the distribution 



of failure times 
the methods wil 



|. If the chords are uncorrelated all 
ive identical results. An advantage 
of approximation (|30|) and its generalization [Eqs. fl2q) 
and (29)] is that Sn and S c i are relatively simple to eval- 
uate; the functions R uu , Rdd and R u d appearing in the 
integral equations Eq. ( [l7|) and ( p3| ) are quite complex 
for a non-zero threshold |p2| . 

To evaluate the approximations for an arbitrary 
threshold we have the following results, 



bi = - + -erf-4=; n c = -e * p 
2 2 V2' 7T 



Sn 



1 f 9{x) dt 



: exp 



(31) 
(32) 
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5cl = y + y erf 



7/9(1 -g) 



g exp 



~i+gj 



erf 



9> 



1 + 9 



(33) 



where 7 = y/-g"{0)and \G\ = ±(1 - g 2 (xj) ~ [g'{x)\ 2 . 
The results for n c |23|| and Sn |13| are well known, and 
the final expression for S c \ can be evaluated using the 
method of Rice as follows. For the level-cut Gaussian 
random field, we have 

S cl = (S(/3 - y(x 1 ))\y'(x 1 )\H((3 - y(x 2 ))) (34) 

The variables w = [wi,w 2 ,w 3 ] = [y(xi), y(x 2 ), y'{xi)} 
have Gaussian distributions with cross correlation ma- 
trix 



gij = (wiWj) => G = 



1 




1 

~9'{x) 





-9'(x) 
-9"(0) 



(35) 



forx = |aj2 — £Ci| andx 2 > Xx- If x 2 < Xi, (y(xi)y' (x 2 )) — 
g'(\x\) but this not effect the final result. To find 5 c i we 
must therefore evaluate 



dwd((3-wx)H(/3-w 2 )\w 3 \ 



(2tt)*|G|! 



(36) 



The integrals extend over all space and the final factor 
is just the joint probability density function of Wi. The 
result is given in Eq. (|3"3|). 

To obtain pi it is necessary to invert Eqs. ( p8| ) and (29). 
This can be done using a short and efficient algo- 
rithm |33| . As previously noted, Pi needs to be known to 
around nine significant figures to achieve four significant 
figure accuracy in the result ]30| ]. To minimize cancella- 
tion errors in the numerators and save one integration, 
we rewrite Eqs. (p^) and (p9f) as 



Pi 



P2 



q/i _ oT 



(37) 
(38) 



where S^i(x) — S c i(x) — \n c . The Laplace transforms 

S'u(s) and 5*Ji(s) on the right-hand side of Eqs. $3% 
and (pq ) can be evaluated using numerical quadrature. 

To check the validity of the independent interval ap- 
proximation, we measure the chord-distribution directly 
from realizations of the thresholded model. This is sim- 
pler (and minimizes finite-size effects) in one dimension. 
A ID random process yi(x) can be obtained from a 
3D GRF 2/3 (r) by taking 7/1(2;) =2/ 3 (r + nx) where r 
is an arbitrary origin and n is a unit-vector with arbi- 
trary orientation. Now yi(x) can be generated indepen- 
dently of 2/3 (r) by using the iD definition for y(x) given 



in Eq. (pq). To ensure that y(x) and y\{x) are statis- 
tically identical they must share g(x). This is true if 
F(k) — 4-7T J k sp(s)ds, where F and p are, respectively, 
the spectral densities of the ID and 3D random fields. 
This shows that F(k) must be a non-decreasing function 
for ID random processes obtained from 3D random fields. 

In the modeling of random media, the following Fourier 
transform pairs [g(x) and p(k) = — F' (fc)/(47rfc)] have 
proved useful 



e -*/C(i + 9) 



4i 



9b = e 
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tt(1 + £ 2 fc 2 ) 2 



/ 7T 



" x/€ (i + ^/0 
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2-jrx/d 
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2tt 2 



tan 1 c_ + tan 1 C4 



c + 



i + ci 



1 + c2 _,; c ± = z(-±k) 



,2tt 
~d 



(39) 

(40) 

(41) 
(42) 

(43) 
(44) 



For a finite number of crossings per unit length (or spe- 
cific surface in three dimensions) it is necessary that 
g'(0) — p3j. For simplicity we restrict attention to the 
following parameters which give 7 = \J—g"{G) = 1/zm -1 : 
(a) £ = 1/im; (b) = \/2pm; and (c) £ = \/2pm, 
d = 47r/v / 6/im. A cross-section of the two-phase medium 
generated in each of the three cases is shown in Fig. |2|(a)- 
(c). We have checked the approximation in the vol- 
ume fraction range <fri € [0.1,0.9]; results for <px=0.2 are 
shown in Fig. ^ (and are typical of those at other vol- 
ume fractions). The independent interval approximation 
is seen to provide remarkably accurate estimates of the 
measured chord-distributions. 

The largest deviations between simulation and the ap- 
proximation are seen for the oscillatory correlation func- 
tion g c (x). We can investigate this "narrow-band" limit 
by taking £ — > 00 in g c {x) which gives, 



9d{x) 



sin 2nx/d 
2-Kx/d 



(45) 
(46) 



The results for p\ at volume fractions 0i=O.2, 0.5 and 0.8 
are shown in Fig. ||, and <Si;i(x) is shown in Fig. || 
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FIG. 2. Chord functions of single-cut Gaussian random 
fields at volume fraction <f>i = 0.2. The symbols are di- 
rectly measured from simulations, and the lines correspond 
to the independent interval approximation. The models [see 
Eqs. (|39[)-(44|)] are shown in the top row with side-length 
40Mm-(a) g a (x) - o; (b) g b (x) - □; (c) g c {x) - A. The 
chord functions correspond to the first-passage time problem 
considered by Rice. 



At 0i=O.5 the approximation is equivalent to that of 
Refs. [P7| , ^8| . The approximation breaks down after one 
wave-length d and actually falls below zero (which is not 
inconsistent with the derivation). This is because the 
process [see Fig. |](a)] has approximately periodic regions, 
extending over several wavelengths, which implies some 
level of correlation between adjacent chords. For exam- 
ple, at </>i=0.5, a chord of length ss is more likely to 
be followed by another of approximately the same length 
than if it were randomly chosen according to the prob- 
ability distribution P2(z). This contradicts the assump- 
tions of the independent interval approximation. The os- 
cillations in the auto-correlation function (Fig. ||) clearly 
reflect appreciable order in the system. 

Note that gd(%) represents the worst-case (or most 
narrow-band) process corresponding to a 3D field since 
Pd(k) oc 6(2?:/ d — k) (ie. an infinitely narrow band pass 
filter). However, the related ID process corresponds to a 
low pass filter (ie. it is not strongly narrow-band). 
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FIG. 3. The phase 1 chord function of a level-cut GRF 
with g(x) = sin(fcr)/(fcr) (k = \/3M m ) a t three different vol- 
ume fractions. The inset shows the microstructure at cf>\ — 0.2 
(side-length 40/im). A ID transect (length 150/^m) of the ran- 
dom field is shown in (a). The horizontal lines corresponds to 
the threshold at each volume fraction. The independent in- 
terval assumption breaks down because the process is nearly 
periodic in some regions. 
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FIG. 4. The normalized two-point correlation function 
[Eq. (^)] of the random model shown in Fig. ^. The strong 
oscillations in S\\(x) correspond to periodic correlations in 
the microstructure. 

This shows that approximations valid for medium and 
wide-band Gaussian random processes are sufficient to 
reproduce the chord functions of 3D models based on 
random fields. 
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The failure of the independent interval approximation 
in this narrow-band limit is not critical for two reasons. 
First, the model g c (x) (with £ > 0, for which the approx- 
imation is reasonable) has been found more relevant to 
physical materials than model gd{x) (e.g., see Sec. 0). 
Second, even in the worst case, the approximation re- 
mains useful out to one wave-length. This may prove 
adequate for material characterization. 



V. EXTENSION TO MORE COMPLEX MODELS 



So far our results have been concerned with the con- 
ventional first-time distributions associated with an arbi- 
trary threshold (/3) of a Gaussian random process. How- 
ever, the single-cut random field model is not sufficiently 
general to model the microstructure of many interesting 
materials. To model the bi-continuous structure of micro- 
emulsions, Berk jll| suggested that phase 1 be defined as 
the region in space where a < y(r) < (3, this is the so- 
called two-cut model [Fig. ||(a)]. The two-cut model has 
also proven useful in interpreting conductivity and perco- 
lation behavior in polymer blends |34}| . Open cell foams 
(e.g. aerogels) and the porous network of sandstones have 
been modeled by the intersection sets of two statistically 
identical (but independent) two-cut fields [Fig. |](b)], and 
closed cell foams may be modeled by the union of two 
such structures [Fig. |](c)] ||. Our method can be sim- 
ply extended to these problems. 

For Berk's Il3[] two-cut model, we have 
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1-.9 2 



xerf 
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(52) 
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1.0 
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FIG. 5. The independent interval approximation compared 
with simulations of the chord functions for three distinct mod- 
els based on level-cut Gaussian random fields: (a) Berk's 
two-cut model - o; (b) the intersection set of two two-cut 
models - □; (c) the union set of two two-cut models - A. 
The side length of the images is 15/^m. 



Using these results we can directly apply the approxi- 
mation for the chord-distributions. We use the field-field 
function gb(x) with Iq — ^/2fj,m and consider a "centered" 
two-cut field (a = —0) at volume fraction 0=0.2. The 
results are shown as circles in Fig. [5] and show very good 
agreement with simulations. 

To evaluate the chord functions of the intersection and 

(49) union sets we first derive their statistical properties. Sup- 
pose f2(x) and Vt'(x) are the indicator functions of two 

(50) independent, but statistically identical, models of ran- 
dom media with properties (f>%, n c , Su and S c %. A new 
model is obtained by forming the intersection set of ft 
and ^ which has indicator function J(x) = fi(x) x ^(x). 
Clearly <j){ = (I) = = eft and n l c = (\I'\) = 

(51) (|fi't& + = + fi|tf'|) = 2<0)(|Q'|) = 2<Mc 



The relation + ft*' 



Ifi'l* + ft|*'| is true ev- 



erywhere expect where the interface of $ and ft in- 
tersect. The contribution of this error to the final re- 
sult is negligible. A similar reasoning can be applied 
to find S(x and S^, as well as corresponding results for 
an analogously defined union set with indicator function 
7(x) = ft(x) + *(x) - ft(x) x *(x). In summary the 
results needed to apply the approximation for the chord- 
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distribution function are 
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(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
(59) 
(60) 



Here the un-superscripted microstructure properties cor- 
respond to the primary models f2 and \& and the super- 
scripted (I or U) functions are to be used in Eqs. fl2S| ) 
and (H). 

Although Eqs. (|53|)-(p0|) are true for any independent 
random models, we restrict attention to the case where 
the primary sets are obtained from Berk's model [see 
Eqs.(p7|)-(50) and Fig. |||(a)]. As above we consider cen- 
tered models (a = — (3) at volume fraction (pi— 0.2 ob- 
tained from random fields with correlation function gb(x) 
(lo = ^/2[im). The results of the independent interval 
approximation are compared with simulations in Fig. |^. 
In general the approximation is excellent. For r < 3/im 
significant deviations (up to 10%) are seen between the 
calculated and simulated values of the chord-distribution 
of phase 2 (p 2 ). 



VI. APPLICATION TO POROUS AND 
COMPOSITE MATERIALS 

To study the properties of a random medium it is im- 
portant to have an accurate model of the microstructure. 
If the physical mechanisms responsible for the evolution 
of the microstructure are not well known (or difficult to 
simulate) an empirically based statistical model may be 
useful §-|l0],[ij||l(|. The level-cut GRF model is well- 
suited to this approach because of its generality: the 
morphology of the model may be "tuned" to some de- 
gree to match that of the random medium. The simplest 
and most common morphological quantities are the den- 
sity (or porosity) and the two-point correlation function, 
which can both be measured from a cross-sectional im- 
age. It is possible to generate a GRF model with approx- 
imately the same statistical properties by an appropriate 
choice of parameters ||,[l2|,[l6| . 

As an example we show a binarized image of a 
tungsten-silver composite |35| along side a single-cut 
GRF model in Fig. The parameters of the model were 
derived in Ref. |36| as follows. The level-cut parameter 
is taken as (3 = —0.84 so that the silver volume frac- 
tion 0i=O.2 is exactly that of the composite [Eq. (31)]. 
The random- field is generated using p c (k) [Eq. (44)]. 
The length scales of the random field £ = 2.15/im and 
d = 13.0^tm are chosen (by a non-linear least squares 




FIG. 6. A binarized image of a silver-tungsten compos- 
ite |3f| (a) compared with a model based on a level-cut 
Gaussian random field (b). The side-length is 99.4/xm. The 
parameters and the model are chosen to reproduce the ex- 
perimental two-point function and chord-distribution func- 
tion (see Fig. Q). 
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FIG. 7. The main graph shows a comparison between 
the chord functions of a silver-tungsten composite (sym- 
bols) and the results of the independent interval approxi- 
mation for a level-cut Gaussian random field model. The 
composite and model are shown in Fig. |^. The inset com- 
pares the experimental and model auto-correlation function 



method) so that the two point correlation function of 
the model matches that of the composite. The theoreti- 
cal and experimental values of Sn(r) (which are practi- 
cally indistinguishable) are shown in the inset of Fig [t]. 
Since the volume-fraction and two-point function do not 
uniquely specify a random microstructure (ie. many dif- 
ferent models may reproduce these morphological quan- 
tities ||), it is necessary to test the results. The chord 
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functions are ideal in this regard as they provide a strong 
signature of microstructure and can be measured from a 
cross-sectional image. The independent interval approx- 
imation and experimental data are compared in Fig. |?]. 
The reasonable agreement between theory and experi- 
ment indicates that the model is capturing important 
features of the tungsten-silver composite. 

A second example is provided by a digitized image 
of Fontainebleau sandstone obtained by X-ray tomog- 
raphy |37|]38| ]. To mimic the granular character of the 
sandstone [Fig. ||(a)] we use a model based on the inter- 
section set of n(=5) single-cut Gaussian random fields. 
The result is shown in Fig. ||(b). To match the porosity 
of the model with that of the sandstone (0=0.154), we 
take P = 0.48 for each of the five primary random field 
models. This corresponds to 4>i — (0.154) 1 / 5 . The ex- 
perimental two-point function is reproduced by choosing 
£ = 51.9/im and d — 272/im in model g c {r) (by a least 
squares method). The independent interval approxima- 
tion for the chord functions is calculated using the re- 
lations 5d = (Sn) n and = n(5ii)"- 1 S' c i [which 
are a straightforward extension of Eqs. fl55|) and (^6|)]. 
The results are shown in Fig. The model is able to 
mimic the the two-point function extremely well, and 
the chord functions with good accuracy. This provides 
evidence that the model is reasonable. 3D images of 
the model and sandstone microstructures are shown in 
Figs. [lO] and O. The sandstone appears more well con- 
nected than the model; the model showing more isolated 
pores. This is actually an artifact of the method used 
to plot the pore-solid interface. An algorithm was used 
to determine that 98.8% of the pore-space in the model 
is connected to the outer faces which compares well with 
99.6% for the sandstone. Therefore the model is also able 
to capture the interconnections of the sandstone pores. 



VII. CONCLUSION 



We have derived a semi-analytic approximation for the 
chord distribution functions (pi and P2) of 3D random 
media. The approximation is based on the assumption 
that successive chord-lengths are uncorrelated. The re- 
sult can be applied to models for which the two-point 
(S11) and ID "surface-void" (S c i) correlation functions 
can be evaluated. The calculation of Su and S c i is gener- 
ally much easier than calculation of p\ and p2- The result 
is exact for Boolean models with convex grains since the 
assumption of independent intervals is true. We have ap- 
plied the approximation to the single level-cut Gaussian 
random field model of random materials. In this case the 
chord functions correspond to Rice's first-passage time 
distribution for random noise. The approximation is very 
accurate for wide-band random fields, but loses accuracy 
after one "wavelength" of the field for extremely nar- 
row band (approximately periodic) fields. Note that a 
narrow-band field corresponds to a low pass (rather than 




FIG. 8. A cross-section of Fontainebleau sandstone (a) 
compared with a model J3(| based on the intersection set of 
five level-cut Gaussian random fields (b). The side-length 
is 2.18mm. The statistical properties of the sandstone and 
model are compared in Fig. 
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FIG. 9. The chord functions measured from a 3D image of 
Fontainebleau sandstone (symbols) compare well with the re- 
sults of the independent interval approximation (main graph) 
for a Gaussian random field model. The inset shows the the 
auto-correlation function 7(2;) = (Su(x) — 0i)/(0i — 4>{). 
Three dimensional realizations of the sandstone and the model 
are shown in Figs. [To] and [ll] respectively. 



a narrow pass) filtered process in ID. 

The result also gives accurate results for Berk's two-cut 
GRF model and other models based on the intersection 
and union sets of level-cut GRF's. This is important for 
generating 3D models of random media using empirical 
information measured from cross-sectional images; the 
two-point correlation function does not necessarily pro- 
vide sufficient information, and good approximations for 
the chord functions are very useful. 
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FIG. 10. 3D representation of Fontainebleau sandstone 
sample obtained by X-ray tomography [ fjT] , ^ . The pore space 
is shown as solid to aid visualization. The side-length of the 
image is 960pim. 

In this context it is possible to apply the approximation 
confidently if the two-point correlation function exhibits 
no (or weak) oscillations. To demonstrate the applica- 
tion of our results we have compared the approximation 
to experimental data obtained from images of a tungsten- 
silver composite and a porous sandstone. 

In order to derive the chord function approximation 
we studied the independent interval process in detail. We 
have shown that the process underlies the derivation, and 
provides useful links between, important results in many 
different fields. The general treatment of the process 
makes clear the relation between various approximations 
for different first-passage times made in signal theory, 
the analysis of component failure and persistence times 
in coarsening. From the expressions for pi [Eqs. ( p8| ) 
and (|29|)1 it is simple to obtain the lineal-path functions 
as Li = £~ 1 {<f)i/s — n c (1 + Pi(sj) /2s 2 ]. For i = 1 this is 
just the "survival probability" in the context of random 




FIG. 11. A sandstone model based on Gaussian random 
fields (cf. Fig. |l^). The independent interval approximation 
for the chord-distribution shows reasonable agreement with 
experimental data (Fig. |^). Many of the apparently isolated 
regions of pore space (shown as solid) are an artifact of the 
plotting procedure. 

processes. Similarly the probability density of "time 
to failure" given that t = falls in a safe region is 
fi = £~ 1 {n c (1 — pi(s)) /2(j>is]. These expressions can 
be inverted in the same way as pt. Due to its apparent 
generality it would be useful to explore the properties 
of the process further. Extensions to include correlation 
between the chord-lengths and the development of a 3D 
analog would be useful future studies. 
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